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^^ ' Abstract: We consider an M/M/oo service system in which an arriving customer is served by the first 

idle server in an infinite sequence Si^ S2, ■ ■ ■ of servers. We determine the first two terms in the asymptotic 

f^ ' expansions of the moments of L as A — > 00 , where L is the index of the server Sl serving a newly arriving 

customer in equilibrium, and A is the ratio of the arrival rate to the service rate. The leading terms of the 
moments show that L/\ tends to a uniform distribution on [0, 1]. 
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1. Introduction 

We consider a stream of customers, with independent exponentially distributed interarrival times, ar- 
riving at rate A to an infinite sequence 51,5*2,... of servers. Each arriving customer engages the server 
Si having the lowest index among currently idle servers, and renders that server busy for an independent 
exponentially distributed service time with mean 1. This stochastic service system, which is conventionally 
denoted M/M/oo, has been extensively studied in the limit A ^ oo; see Newell [N]. We shall be interested in 
a question mentioned only tangentially by Newell: what is the distribution of the random variable L defined 
as the index of the server Sl serving a newly arriving customer when the system is in equilibrium? Newell 
[N, p. 9] states that L "is approximately uniformly distributed over the interval" [1, A], basing this assertion 
on the approximation 

Pr[L>/]« < ^" A' if'<-^' (1.1) 

i 0, if / > A. 

But no error bounds are given for this or other approximations stated by Newell, and not even the fact that 
the first moment has the asymptotic behavior 

Ex[L] ^ ^ (1.2) 

that it would have under the uniform distribution is established rigorously. Our goal in this paper is to give 
a rigorous version of (1.1) that will suffice to establish not only (1.2), but also the next term, 

Ex[i] = ^ + ilogA + 0(l), (1.3) 

and more generally 

Ex[L-]^^::^ + -^"'^^"g^+0(A-^) (1.4) 

for 771 > 1. In particular, we have 

Var[L] = Ex[L2] - ^x[Lf 

Since the interval [0, 1] is bounded, formula (1.4) shows that the m-th moment of L/X tends to l/{m + 1) as 
A — > oo for all TO > 1, and thus suffices to show that the distribution of L/X tends to the uniform distribution 
on the interval [0, 1]. We note that a problem that is in a sense dual to ours (finding the largest index of 
a busy server, rather than the smallest index of an idle server) has been treated by Coffman, Kadota and 
Shepp [C]. 

The key to our results is the probability Pr[i > /], which is simply the probability that the first I servers 
5i , . . . , 5; are all busy. It is well known that this probability is given by the Erlang loss formula 

Pr[L > f - ''^'■ 
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D, 



where 



0<k<l ^ ' 



(see for example Newell [N, p. 3]). The sum Di can be expressed as an integral, 

(see for example Newell [N, p. 7]), and most of Newell's analysis is based on such a representation. But we 
shall work directly with the expression of Di as the sum in (1.5). 

We shall divide the range of summation in (1.5) into two parts. The first, which we shall call the "body" 
of the distribution, will he < k < Iq — X — s, where s — ^/X. The second, which we shall call the "tail", 
will be I > lo- In Section 2, we shall derive an estimate for Pr[L > I] in the body, and in Section 3, we shall 
derive an estimate for the tail. In Section 4, we shall combine these estimates to establish (1.4). 

2. The Body 

In this section we shall establish the estimate 

Pr[L > /] . (1 - l/X) + 3^^^ + O (i) + O {y^^^) (2.1) 

for I < Iq = X — s, where s — Va. We begin by using the principle of inclusion-exclusion to derive bounds 
on the denominator Di. 

We begin with a lower bound. Since 



fc-i 




we have 







D V '^^-')' 


■ ■(i-k + i) 


0<fc<( 


Xk 




1 ^ /k\ fl^'^ 

X ^ \2j \Xj 


y ^-V 


i + o{ii/xy) 

i-i/x ■ 



For the first sum we have 



We note that the logarithm of (l/XY has a non-negative second derivative for I > 1. Thus (l/XY assumes 
its maximum in the interval < ^ < ^o for I — 0, I = 1 or I = Iq. Its values there are 0, 1/A and 
(1 - s/X)^^" = (1 - 1/VX)^^^ < e"^+\ respectively. As A ^ oo, the largest of these values is 1/A, so 
wc have 0{(l/Xy) = 0(1/A) for < ^ < Iq- Thus the first sum is 

y- fl\''_ 1 + 0(1/A) 



For the second sum we have 



XJ l-l/X 



V /^^ ^'^''"'- i+o(/2(VA)') 



,2/ \X (l-^/A)3 



The logarithm of /^(//A)' has a non-negative second derivative for ^ > 3, so an argument similar to that used 
for the first sum shows that 0{f{l/Xy) = 0(1/A) for < ^ < lo- Thus we have 

oh V2^ V a; (1 - //A)3 

and the lower bound 

1 + 0(1/A) 1 + 0(1/A) 

^'- 1-i/x ~x{i-i/xr- (^-^^ 

For an upper bound, we have 



\0<j<k-l J \o<i<j<k-l 
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(because Eo<^<,<fe-i ».? = ^Eo<,<fe-i.?j - Eo<,<fe-i j'j /2 < (Eo<,<fc-i.?j /2 = (2) /2)- Thus 
we have 



,A/ A ^^ V2/ V A/ 2A2 ^^ V2 



For the third sum we have 

e(:)(t) ^e 






,2/ VA/ -^^V2/ VA 

0<fc<( ^ ^ ^ ^ fe>0 



and thus the upper bound 

^^ ^ 1 + 0(1/A) _ 1 + 0(1/A) ^^f 1 



1-l/A A(l-VA)-3 VA2(1-VA)5 

Combining this upper bound with the lower bound (2.2) yields 

^ _ 1 + 0(1/A) 1 + 0(1/A) ^ ^^ 1 



1-VA A(1-VA)3 VA2(1-VA)-^ 

To obtain Pr[L > Z], we take the reciprocal of Df. 



1-Z/A A(l-//A)3 VA2(1-VA)^ 

(l + 0(l/A))(l-VA)(l-3^(^^ + 0( 



A(l-//A)2 VA2(1-VA)^ 

- (1 + 0(1/A)) (1 - VA) (1 + 3^(^^ + O {-^^^^^ 

-i' + ^^y^^) ((i-^/^) + A(r^ + KA^(r^))- 

Observing that 0(1/A) (1 - l/X) = 0(1/A) and 0(1/A)/A(1 - l/X) = 0(1/A2(1 - l/Xf), we obtain (2.1). 
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3. The Tail 

In this section we shall establish the estimate 



PT[L>l]^0{e-^X'/U) (3.1) 

for I > X — s, where s == vA. To obtain an upper bound on Pr[L > I], we obtain a lower bound on D;. We 
have 

'" ^ (i-ky.xi^ 

0<k<l ^ ' 

- [X-s\l A'- L^-^J ^ ■ " ^ [A - 2sJ ! A'- L^-2-J ' ^^'^^ 

because / — [A — ,sj > ^ — (A — s) > by assumption and [A — 2sJ > for all sufficiently large A. There 

are [A — 2s\ — [A — 2s\ + 1 > s terms in the sum (3.2). Furthermore, the smallest of these terms is the 

last, because its denominator contains factors of A where the preceding terms contain factors smaller than 

A. Thus we have 

sl\ 
Di> 



[A-2sJ!A'-L-^-2-J 



^ — n ^n 



for all n > 1 (because the trapezoidal rule underestimates the integral J^ log x dx of the concave function 
logs). This estimate yields 



We have 



Di > , . (3.3) 

" e ^[A-2sJ [A - 2sJ l>-'^^\ A'- L^-2''J 



LA-2sjL^-2.J < (A-2s)L^-2sj 

= AL^-2sJ(^_2s/A)L^-2.J 

< AL^-2sJ (^ _ 2s/A)^-2.-i 

< AL^^^^J g(-2s/A)(A-2s-l) 
^ ;)^L^-2sJ g-2s+4sVA+l 



and 



Substituting these bounds into (3.3) yields 



\/lX-2s\ < s. 
Di> 



e^A'' 
Taking the reciprocal of this bound yields (3.1). 



4. The Moments 

In this section we shall use (2.1) and (3.1) to prove (1.4). We write 

A^ii) ^ r - {I ~ ir 

for the backward differences of the ?Ti-th powers of I. Then partial summation yields 

= ^A,„(OPr[L>/] 
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J2 m r-i Pr[L > /] + O XI ^""^ P^[^ > '] ("^-1) 



i>o \ i>0 



This formula shows that we should evaluate sums of the form 

T„ = Xr Pr[L>;]. (4.2) 

We shall show that 

(n+l)(n + 2) 2 

Substitution of this formula into (4.1) will then yield (1.4). 

We shall break the range of summation in (4.2) at Iq — \ — s, where s = a/A, using (2.1) for < ^ < Zq 
and (3.1) for I > Iq. Summing the first term in (2.1), we have 



E ^"(l-VA) = y E (A;"-?"+^) 



A 
o<;<io o<(<io 



1 ^^A(A«+i-(n + l)A"5) , ^,^„,^ ^A"+2-(n + 2)A"+i,s , ^,^„+,, 



0(A") - ^^-^ + 0(A'^ 



AVV "^+1 /V n + 2 

\n+l 

0(A"). 



(n+l)(n + 2) 

Summing the second term in (2.1), we have 

,r- J^_ V- (A-fc)" 
^ X-l^ ^ k 

0<l<la s<k<X 

s<fc<A ^ 

= A"log- + 0(A") 

s 



where we have used X]i<A;<ri ^/^ ^ logn + 0{1). Summing the third term in (2.1) of course yields 0{X" 
Summing the last term in (2.1), we have 

A V ^" = A V (^~^)" 
^ (A-/)3 ^ P 

0<l<lo ^ ' s<k<\ 

< A"« S ^ 

s<fe<A 



fc3 

k>s 



= 0(A"), 
where we have used X]fc>n 1/^^ = 2/n^ + 0(l/n3). Combining these estimates, we obtain 

Finally, summing (3.1) we have 

= 0(A"), 
because the summation on the right-hand side is the n-th moment of a Poisson random variable with mean 
A, which is a polynomial of degree n in A. Thus 



^r Pr[L>/] =0(A") 



Combining this estimate with (4.4) yields (4.3) and completes the proof of (1.4). 

5. Conclusion 

We have obtained the first two terms in the asymptotic expansions of the moments of _L as A — > oo . An 
obvious open question is whether one can obtain a complete asymptotic expansion, or even just the constant 
term in (1.3) and the corresponding terms in (1.4). While our estimates for the contributions to the 0{\) 
term in (1.3) could be improved (for example, by a better choice of the parameter s), it is clear that new 
techniques will be needed to obtain an error term tending to zero. 
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